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Nonlinear conductivity effects are studied experimentally and theoretically for thin samples of 
disordered ionic conductors. Following previous work in this field the experimental nonlinear con- 
ductivity of sodium ion conducting glasses is analyzed in terms of apparent hopping distances. Values 
up to 43 A are obtained. Due to higher-order harmonic current density detection, any undesired 
effects arising from Joule heating can be excluded. Additionally, the influence of temperature and 
sample thickness on the nonlinearity is explored. From the theoretical side the nonlinear conductiv- 
ity in a disordered hopping model is analyzed numerically. For the ID case the nonlinearity can be 
even handled analytically. Surprisingly, for this model the apparent hopping distance scales with 
the system size. This result shows that in general the nonlinear conductivity cannot be interpreted 
in terms of apparent hopping distances. Possible extensions of the model are discussed. 

PACS numbers: 66.10.Ed, 66.30.Hs, 61.43.Fs 
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I. INTRODUCTION 

One important method for improving the properties 
of solid ionic conductors with regard to applications in 
microbatteries, fuel cells and electrochromic devices is 
the preparation of thin films [JEHSHHSH. Both 
the electrical resistance and the mechanical stiffness of 
the film decrease with decreasing thickness. The decrease 
of the electrical resistance is important for improving the 
power density of microbatteries, the efficiency of fuel cells 
and the switching time of electrochromic devices. A low 
mechanical stiffness leads to a better processibility of the 
film. 

For the integration of a thin film into a lithium micro- 
battery, an electrochemical window in which the film is 
stable from 0-5 V is required Q, . The application of 
5 V to samples with thicknesses of 100 nm and less leads 
to high electric field strengths in the samples, which re- 
sult in field-dependent ion transport effects. Generally, 
the ionic conductivity of solid electrolytes increases with 
increasing field strength H [H EJ III Thus, field- 
dependent ion transport is of potential interest for im- 
proving the applicability of thin-film electrolytes. 

Mathematically, the field-dependent electrical proper- 
ties of thin solid electrolyte samples can be described by: 



turns out that 



jdc{Edc) = Cl,dc Edc + &3,dt 



E dc 



£3 



(1) 



Here, the linear dependence of the dc current density jdc 
on the dc electric field Ed c at low field strengths is char- 
acterised by the low-field dc conductivity a\^ci while the 
nonlinear dependence at high field strengths is charac- 
terised by the higher-order conductivity coefficients oz,dc-> 
o~5,dc etc. 

A first step to interpret the nonlinearity of the con- 
ductivity is to consider a simple regular hopping model 
with distance a between adjacent sites. For this model it 



jdc(Edc) oc smh(q(3aE dc /2) 



(2) 



with P = l/fc^T. Thus, in the framework of this model, 
it is possible to extract the hopping distance a from field- 
dependent electrical data. Although the model is too 
simple to provide a realistic description of ion conduc- 
tion in disordered materials, such as glasses and polymer 
electrolytes, the jdc(Edc) curves of many real ion conduc- 
tors can be reasonably well fitted by Eq. ©. However, 
the values obtained for the hopping distance a are much 
larger than typical distances between neighboring sites 
in ionic conductors. For instance, in the case of ion con- 
ducting glasses, values in the range from 15 A to 30 A 
have been found 0, E3, EH 

These experimental results are based on measurements 
using dc electric field. A major drawback of this method 
is the lack of information about Joule heating effects. 
Joule heating may lead to an increase of the sample tem- 
perature, resulting in an increase of the ionic conduc- 
tivity. In contrast, the application of ac electric fields 
allows for a direct differentiation between nonlinear ion 
transport and Joule heating. Nonlinear ion transport 
leads to higher harmonic contributions to the current 
density spectrum, while this is not the case for Joule 
heating [lj|. The field dependence of the higher har- 
monic contributions can be used to determine values for 
the higher order-conductivity coefficients 03^0 cr5.de etc. 
|14j . Thus, by means of nonlinear ac conductivity spec- 
troscopy it is possible to differentiate between the higher- 
order conductivity coefficients, while the application of 
dc fields yields only one jdc(Edc) curve. 

The values a\^c and a^^c can be used to define an 
apparent hopping distance: 



app 



= 2^ 4c /{a 14c q 2 /3 2 ) 



(3) 



This definition of a app implies that for the regular hop- 



ping model, the simple relation a app = a holds. However, 
two of the present authors have shown that for differ- 
ent sodium ion conducting glasses, the apparent hopping 
distances a app are in a range from 39A to 55A Hi. 
Thus, the calculation of a app by means of Eq. leads 
to higher values than fits of jdc(Edc) data by means of 
Eq. Q. This implies that Eq. J5J does not provide an 
exact description of the nonlinear conductivity of ionic 
conductors. This is confirmed by the observation of neg- 
ative values for a 5 ^ c whereas the validity of Eq. @ 
implies positive values for Osse- 
in this paper, we rationalize these experimental results 
by considering simple hopping models with strong site 
and barrier disorder. We find that in the framework of 
such models, the apparent hopping distance a app is, in- 
deed, considerably larger than the distance a between 
adjacent sites. However, a app shows an unexpected de- 
pendence on the thickness of the model systems in the 
direction of the applied electric field. In order to check 
whether also real ionic conductors exhibit such a thick- 
ness dependence of their nonlinear electrical properties, 
we have carried out nonlinear ac conductivity measure- 
ments on ion conducting glass samples of different thick- 
ness. We compare theoretical and experimental results, 
and we discuss implications of these results for further 
theoretical and experimental work. 

II. THEORY 

First, we consider a general ID-hopping model with 
site and barrier disorder; see Fig. ^ The individual N 
sites are characterized by energies Ei. With an external 
field the hopping rates for leaving site i are r^± where 
the sign denotes whether the hopping is along (+) or 
opposite (-) to the field. The field dependence of r^-t 
can be expressed as Ti t ±(u) = exp(±u)7i j ± with 

u = qf3aE dc /2. (4) 

Correspondingly, 7^-1- are the hopping rates without ap- 
plied field. Detailed balance requires 7i i _i_/7i+i,- = 
exp(— (3(Ei + i — Ei)). To get a stationary long-time so- 
lution for the current, periodic boundary conditions are 
used. Physically, this reflects the interaction of both elec- 
trodes via the voltage source. 

The time evolution of the single-particle problem is 
governed by the rate equations 

(d/dt) Pt (t) = -r it+ (u) Pi (t) - Ti^(u)pi(t) 

+T i _ li+ (u)p i -i(t) + r 4+ i _(u)p i+ i(i)(5) 

The stationary long-time solutions are denoted as p°° (u) . 
The argument expresses their dependence on the electric 
field. Finally the current jdc(u) in the long-time limit 
can be written as 

3dc{u)/{qa) =j(u) =pf(u)r i)+ (u) -p^. 1 (u)T i+1 _(it). 

(6) 



2 




1 1 1 1 1 > 

a (N-2)a (N-1)a x 

Figure 1: Sketch of the ID hopping model, used in the theo- 
retical analysis. Periodic boundary conditions are employed. 
The model is point symmetric. 

Note that the current between all pairs of sites is iden- 
tical. Therefore j(u) and thus jdc(u) do not depend on 
the index i. 

This property can be used to solve the system of equa- 
tions analytically for arbitrary N; see Refs. [T^, [l(| H3, 
Hs| for detailed information. 

Using a notation which is appropriate for our further 
analysis the solution can be written as 

j(u) = (exp(iVu) - exp(-Nu))/B (7) 

with 

N-i 

B =Y^ A m exp((l - N + 2m)u). (8) 

m=0 

The disorder is contained in 

N—i 

A m =Y / f( m ' i ) (9) 

i=0 

with 

f(m,i) = exp(-/3(£ i+m - ^))/ 7l ,_. (10) 

In case that the disorder of adjacent sites is uncor- 
rected one obtains (^4i) = ... = (Ajv-2) = b\ and 
(A ) = (Ajf-i) = b , where the brackets denote the dis- 
order average. Please note that (A m ) as a sum of TV 
terms trivially scales like N . 

In what follows we consider the situation that the en- 
ergy landscape is point symmetric. In this way one can 
guarantee that in analogy to the experimental situation 
the current only depends on odd powers of the electric 
field. In particular one strictly has A m = A^^i^ m . 

Here we are interested in a Taylor-expansion of j(u) 
with respect to u. Expansion of the numerator and the 
denominator yields 

- _ 2Nu + N 3 u 3 /3 + ... 

J[U) E m A m + (1/2) £ m A m (-N + 1 + 2m)2 u 2 + ... ■ 

(11) 



Reordering of terms yields up to order u 2 
2Nu 



l + u 2, NV6 _E m A m (-N + l + 2m) 



2\ 1 



(12) 



From this the low-field conductivity a\^ c can be written 
as: 



01, dc 

q 2 /3a 2 



N 



Em ^ 



(13) 



The brackets denote the average over the disorder. 

In a first step we use the property that in the thermo- 
dynamic limit iV — ► oo one has 



1 



Em Ar, 



1 



(Em Am) 



(14) 



This relation is derived in the Appendix. 

With the same arguments one can replace in the dis- 
order average of j(u) in Eq. (fT2"|) all terms ^ m Am by 
(£ m i ra ). Combining Eq. J3J) with Eq. ijT2|) and using 
this preaveraging of ^2 m A m one obtains in the limit of 
large N the following relation for a 2 : 



2 I 2 
a appl a 



N 



(Em A-m) 

Straightforward calculation yields 

(J2An) = Nb 1 +2(b Q -b 1 ) 



(15) 



(16) 



and 

(J2A m (-N + l + 2m) 2 ) = (1/3M7V 3 - N) 
m 

+ 2(b - h)(N - 1) 2 (17) 

Expansion of the denominator with respect to 1/N finally 
yields 



alja 2 = 1 + 4JV(6i - b )/h 



(18) 



This shows that for large N one finds the scaling relation 



Kpp oc N . 



What is the origin of the ./V-dependence? In the case 
of vanishing disorder one has pf^iu) = 1/N for all u. 
Thus the M-dependence of the current j(u), see Eq. (JBJ, 
exclusively results from the trivial u-dependence of the 
rates I\ ,±(u) and gives rise to a 2 pp /a 2 = 1. In contrast, 
for disorder one also obtains contributions from p°°(u). 
After solving the N linear equations, characterizing the 
stationary case, the expansion of p°° (u) with respect to u 
contains terms proportional to N, which are responsible 
for the iV-dependence of a 2 app /a 2 . For the N-dependence 
of the <Jj,dc one finds (J\,& c oc TV -1 and o-&,dc oc N , thus 
giving rise to a\ pp oc cr 3 ,d c / '^i,dc oc N. The scaling of 
cti^c simply means that the conductivity is proportional 
to the particle concentration, i.e. 1/N. Actually, it is 
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Figure 2: Comparison of the numerical solution for a apP with 
the asymptotic analytical solution. 



also possible to calculate j(u) for fixed u in the thermo- 
dynamic limit, thereby revealing non-analytical behavior 
for u = 0. This discussion, however, is beyond the scope 
of the present paper. 

It may be instructive to calculate bo and b\ for different 
types of disorder. In case of a random barrier model 
for which all Ei are constant one obtains bo — b\ and 
thus a 2 pp /a 2 = 1. In case of pure energetic disorder we 
choose a model for which Ei is either e or 0, both with 
50% probability. Then one finds numerically in the low- 
temperature limit that bo ~ 0.24exp(/?e) and b\ « 
0.37exp(/3e). For large N one thus gets a 2 exp /a 2 « 1.5N. 

In Fig. we show the asymptotic result Eq. 1181) 
for a 2 pp /a 2 together with the exact result, obtained from 
averaging the current in Eq. l|12fl and finally calculating 



app 



/a 2 from Eq. One can see that for N > 40 the 
large-N limit works very well. For smaller N the appar- 
ent jump length is smaller than expected. We mention 
in passing that all data points have been obtained from 
averaging over 10000 realisations of the disorder. This 
large number is necessary because numerically it turns 
out that the standard deviation for the distribution of 
a 2 pp /a 2 is of the same order as a 2 pp /a 2 itself (for all N). 
The main conclusion is that for larger values of TV also 
large values of a 2 pp /a 2 can be realized. Furthermore we 
have included data for u = 2 in Fig. [21 Decrease of 
u corresponds to an increase of temperature. In agree- 
ment with the experimental data (see below) one finds 
that a 2 pp /a 2 decreases with increasing temperature. Of 
course, for u = 0, one obtains a 2 pp /a 2 — 1. Actually, as 
seen from Fig. [21 the asymptotic result is not strictly lin- 
ear as naively expected from Eq. i|18fl . The reason is that 
for small N also the values of bo/N and bi/N somewhat 
depend on TV due to finite-size effects. 

So far we have considered a single particle in a dis- 
ordered potential. In recent work it has been shown 
HUMES that it is a more realistic picture to consider 
the dynamics of a vacancy rather than a particle. Not 
surprisingly, the new set of equations is very similar and 
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Figure 3: Dependence of a apP on the size parallel and orthog 
onal to the electric field. 



Figure 4: Apparent jump distance of mobile Na + ions, a app , 
in NCAS25 glass samples with different thicknesses. 



the results for a 2 p P /a 2 turn out to be identical. 

A more serious deviation from the experimental sit- 
uation is the dimensionality. To which degree can the 
results for the ID model be generalized to systems with 
more than one dimension? To clarify the influence of 
the dimension we have extended the model to two di- 
mensions, also using periodic boundary conditions in the 
second dimension. We define N Q as the number of sites 
orthogonal to the field. The set of rate equations for the 
ID case can be easily generalized. Solving for the sta- 
tionary long-time solution is equivalent to solve a system 
of linear equations with N ■ N Q variables. Unfortunately, 
it is not possible to solve the multidimensional case ana- 
lytically in the way it has been done for the ID system. 
Therefore the solution is purely numerical. After aver- 
aging over the disorder one obtains the results, shown in 
Fig. [3J Two important conclusions can be drawn from 
these results: (i) The scaling a 2 pp /a 2 oc N also holds for 
two-dimensional systems, albeit only for somewhat larger 
N (as can be seen for 7V = 20). (ii) The absolute value 
of ciapp/a 2 strongly decreases with increasing N a . Com- 
paring N Q = 10 and N a = 20 one may speculate that 
a 2 a pp/a 2 cx N/N Q . 



III. EXPERIMENTAL 

For the nonlinear ac conductivity measurements, we 
prepared samples of the glasses 0.25 Na 2 • 0.096 CaO • 
0.062 A1 2 3 • 0.592 Si0 2 (NCAS25) and 0.25 Na 2 • 0.75 
Si0 2 (NS25) with different thicknesses. The preparation 
of the NCAS25 glass is described in Ref. 44]. The NS25 
glass was prepared by using the same procedure. 

The experimental setup for nonlinear conductivity 
spectroscopy and the high-voltage measurement system 
are described in detail in Ref. |l4|. An important as- 
pect of our method is the utilisation of nonblocking elec- 
trodes consisting of highly conducting NaCl solutions. 
Thereby, we avoid (i) electrochemical reactions at the 



electrode/sample interfaces and (ii) electron or hole in- 
jection into the sample. The measurements were carried 
out in a frequency range from 10 mHz to 10 kHz with a 
maximum voltage amplitude of 500 V. 



IV. RESULTS 

For the NCAS25 glass, the third-order conductivity 
coefficients were derived from an analysis of the 

higher harmonic components of the current density spec- 
trum as described in Ref. [l4|- At low frequencies, both 
the low-field conductivity a\ {v) and the third-order con- 
ductivity coefficient 03(1/) exhibit plateaus, the plateau 
values being identical to the bulk dc values, cri,dc and 
03,<2c 44|. From these dc values, the apparent jump 
distance a app were calculated using Eq. J3J). A careful 
analysis of the frequency dependence of the conductivity 
data revealed, however, that in the low-frequency region, 
sample/electrode interface impedances have a weak in- 
fluence on the conductivity. This will be discussed later 
in more detail for the NS25 glass where the interfacial 
impedance effects are more pronounced. For the deter- 
minination of a\^ c and a-^ dc of the NCAS25 glass, the 
analysis was limited to a frequency regime where the 
interfacial impedance was less than 1% of the overall 
impedance. Due to the narrower frequency regime, the 
values obtained for a app were a few percent lower than 
those published in Ref. ^4|. In Fig. 0|we show a plot of 
CLapp versus temperature T for two NCAS25 glass samples 
with thicknesses 65 /im and 85 /im. Within the exper- 
imental error, we do not detect a significant thickness 
dependence of a app - 

In the case of the NS25 glass, we carried out measure- 
ments on samples with thicknesses 83 ^m, 100 fim and 
120 /im. Unfortunately, for all samples, the higher har- 
monic current density did not exhibit well-defined 
low-frequency plateaus. As an example, we show in Fig. 
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Figure 5: Base current density j'(u) (closed symbols) and 
higher harmonic current density j'(3v) (open symbols), mea- 
sured after applying an ac field with amplitude Eq — 4f.6 
kV/cm to the NS25 glass. Both Fourier components of the 
current density are normalised by the field amplitude Eo and 
are plotted versus frequency v at two different temperatures, 
T = 303 K and 3f3 K. In addition, the quantity j (3v)/Eo is 
multiplied by a factor of -4. 



results for the base current density j'(y) (closed sym- 
bols) and for the higher harmonic current density j'(3v), 
measured after applying an ac field with amplitude Eq = 
41.6 kV/cm to a sample with thickness d = 120 /Ltm. Both 
j'(y) and j'(3u) are normalised by the field amplitude Eq 
and are plotted versus frequency at two different temper- 
atures T = 303 K and 313 K. While a bulk dc plateau in 
j'(y) is clearly detectable, the higher harmonic current 
density j'(3v) exhibits a frequency dependence over the 
entire frequency range. Therefore, it was not possible to 
obtain accurate values for f73.d c of this glass. 

The reason for the absence of wcll-dcfindcd dc plateaus 
in j'(3u) are most likely electrical impedances at the sam- 
ple/electrode interfaces. Although we use highly con- 
ducting liquid NaCl solutions as electrodes, we detect sig- 
nificant interfacial impedances at low frequencies, which 
are related to additional barriers for the diffusion of the 
sodium ions from the glass samples into the NaCl solu- 
tion and vice versa. As an example, we show in Fig. 
plots of the imaginary part of the impedance and of the 
real part of permittivity versus frequency for the NS25 
glass. The data were obtained at T = 313 K and differ- 
ent applied voltages. At high frequencies, the impedance 
and permittivity data are determined by ionic conduc- 
tion in the bulk and exhibit only a weak voltage depen- 
dence. However, below 1 Hz, additional impedance con- 
tributions from the sample/electrode interfaces are de- 
tected. These interfacial impedances decrease strongly 
with increasing voltage. However, even at the highest 
applied voltages, interfacial impedance contributions are 
still present. 
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Figure 6: Frequency dependence of the imaginary part of 
the impedance, Z" and of the real part of the permittivity, 
e', for the NS25 glass. The data were taken at T = 313 K 
and different voltage amplitudes. Above 1 Hz, the spectra are 
determined by bulk ionic conduction, while below one 1 Hz, 
sample/electrode interface impedances govern the frequency 
and field dependence of the electrical data. The interfacial 
impedances are related to additional barriers for the diffu- 
sion of the sodium ions from the glass samples into the NaCl 
solution and vice versa. 



DISCUSSION 



Although we did not detect a significant thickness de- 
pendence of the apparent hopping distance a app for the 
NCAS25 glass, our experimental results do not com- 
pletely rule out such a thickness dependence. For an 
unambiguous prove or disprove we will have to study 
samples with a broader range of thicknesses. In the 
case of the NS25 glass, the thickness was varied by 
about 50%, however unfortunately, accurate a^ dc val- 
ues for this glass could not be obtained, most likely due 
to sample/electrode interface impedance contributions to 
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It is remarkable that the interfacial impedances persist 
up to voltages of several hundred volts. A voltage drop 
of the order of 100 V in an extremely thin interface re- 
gion should lead to enormously high electric fields which 
should pull the ions over the interfacial barriers. When 
we assume that the thickness of the interfacial region is 
of the order of 10 nm, a voltage drop of the order of 100 
V in this region leads to field strengths of the order of 100 
MV/cm. This is even far above the electrical breakdown 
strengths of ion conducting glasses. 

From this we conclude that the interfacial regions 
causing the interfacial impedance effects must be much 
thicker than 10 nm. A possible explanation is the build- 
up of a resistive surface layer due to water corrosion. It is 
thinkable that water diffuses from the NaCl solution into 
the glass and and changes the chemical structure close to 
the surfac e by breaking Si-O-Si bonds and forming Si-O- 
H groups [22| • The thickness of such chemically modified 
surface layers may be much larger than 10 nm [22j. This 
assumption would explain the more pronounced interfa- 
cial impedance effects in the case of the NS25 glass. It is 
well known that the chemical corrosion of simple alkali 
silicate glasses is much faster than that of technical sil- 
icate glasses containing additional alkaline-earth oxides. 
If the assumption is correct, the interfacial impedances 
should be avoidable by using non-aquaous salt solutions, 
based for instance on glycerol. This will be the subject 
of further experiments. 

From a theoretical point of view, the observation of 
large values of a app has been reproduced by a simple 
model, describing the hopping dynamics of a particle 
(or vacancy) in a disordered model potential. In the 
ID case it is even possible to solve the problem ana- 
lytically. This simple and well-studied model contains 
the surprising property that a app scales with the thick- 
ness. Identification of a app = (24cr 3i( ; c /g 2 /^oi^c) 1 / 2 as 
some kind of effective hopping distance was motivated by 
the solution of the regular model. The scaling property 
a app oc aN 1 / 2 for the disordered hopping model directly 
implies, however, that this interpretation is not justified. 
Formally, this can be traced back to the fact that the 
field-dependence of the current in the disordered case is 
not only goverened by the new transition rates T^± but 
also by the modified populations; see Eq. Rather for 
this model the value a app /(aN 1 ^ 2 ) expresses information 
about the nature of the disorder; see Eq. ijTHjl. Thus, the 
experimental observation of a app « 50A does not con- 
tradict the generally accepted idea that typical hopping 
distances in ion conductors should be smaller than, let's 
say, 10 A. 

VI. CONCLUSIONS 

Combining a detailed experimental analysis with a 
thorough model calculation several new aspects have 



been elucidated, (i) The large values for a app in thin 
samples is not related to Joule heating, but is indeed a 
reproducible property of the ion dynamics, (ii) These val- 
ues increase with decreasing temperature, (iii) The large 
values of a app as well as the temperature dependence can 
be rationalized by a simple hopping model, (iv) In gen- 
eral, a app cannot be interpreted as an apparent hopping 
distance. 

Many new and important questions have emerged from 
the present results. (1) What is the real physical inter- 
pretation of o app ? In particular it would be helpful to 
relate this quantity to equilibrium properties of the sys- 
tem |23|,|^J. (2) Is is possible to refine the model to get a 
limiting value for a app ? One may speculate that consid- 
eration of interaction effects between different particles 
or holes, respectively, may yield such an upper limit. (3) 
Do large electric fields also modify the properties of the 
network and thus indirectly also the ion dynamics? The 
latter effect has not been included in the modelling. In 
principle this question could be answered by performing 
appropriate molecular dynamics simulations of ion con- 
ductors. (4) Are experimental values for a app different for 
crystalline materials as compared to disordered materi- 
als? The hopping model, analyzed in this work, might 
suggest some difference. 

Appendix 

We consider the covariance of A^ and A m which can 
be written as 

cov(A k ,A m ) = ^[(/(fc,i)/(m,i)) - (/(M)X/(m, j)>]. 
%,i 

(19) 

Closer inspection of the f(k,i) (see Eq. (|1U|> ) shows that 
for extreme disorder the first term is particularly large if 
i = j. This can be seen from the fact that without the ad- 
ditional exponential factor in Eq. HlOfl one would strictly 
ha,ve(f(k,i)f(m,j))-{f(k,i)){f(m,j)) = Ofori ^ j. As 
a consequence one can write cov(Ak, A m ) as a sum of ./V 
terms of similar magnitude such that cov^Ak, A m ) oc N. 
Thus the variance of ^2 m A m which can be written as 
Sfc m cov (Ak, A m ) scales like N 3 . This has to be com- 
pared with the scaling Em^™) 2 x N 4 - Thus for large 
N the standard deviation of A m is much smaller than 
the value of J2 m A m itself. This justifies Eq. (|HJl . 
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